1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 300,000 iteration;
- 200,000 burn-in;
- 100 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.07
## LDevsum 1.45 2.35
## dh0 1.09 1.31
## dh1 1.09 1.32
## dl0 7.54 16.72
## dl1 1.30 1.93
## dl2 6.79 15.95
## dl3 3.69 7.88
##
## Multivariate psrf
##
## 4.97
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1216.28632 | 22069.538 |
| DIC3 | 1156.95606 | 22120.978 |
| PWAIC | 42.68576 | 234.347 |
| WAIC | 1183.68816 | 22141.168 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 300,000 iteration;
- 200,000 burn-in;
- 100 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 0.999 1.00
## LDevsum 1.001 1.01
## dh0 1.022 1.10
## dh1 1.028 1.12
## dl0 1.651 3.04
## dl1 1.401 2.27
## dl2 1.257 1.95
## dl3 1.446 2.33
##
## Multivariate psrf
##
## 1.47
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1213.74695 | 22043.9831 |
| DIC3 | 1157.49543 | 22121.3879 |
| PWAIC | 43.02139 | 235.2914 |
| WAIC | 1184.60205 | 22142.1846 |